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Abstract 



We present a unified description of nuclear deep inelastic scattering (DIS) 
over the whole region < x < 1 of the Bjorken variable. Our approach is 
based on a relativistically covariant formalism which uses analytical properties 
of quark correlators. In the laboratory frame it naturally incorporates two 
mechanisms of DIS: (I) scattering from quarks and antiquarks in the target 
and (II) production of quark-antiquark pairs followed by interactions with 
the target. We first calculate structure functions of the free nucleon and 
develop a model for the quark spectral functions. We show that mechanism 
(II) is responsible for the sea quark content of the nucleon while mechanism 
(I) governs the valence part of the nucleon structure functions. We find that 
the coherent interaction of qq pairs with nucleons in the nucleus leads to 
shadowing at small x and discuss this effect in detail. In the large x region 
DIS takes place mainly on a single nucleon. There we focus on the derivation 
of the convolution model. We point out that the off-shell properties of the 
bound nucleon structure function give rise to sizable nuclear effects. 

PACS numbers: 13.60.Hb, 25.30.-c 
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I. INTRODUCTION 



Deep-inelastic lepton scattering (DIS) on nuclei is a powerful tool to investigate the 
quark-gluon structure of nucleons in a nuclear environment. The cross-section of this 
process is studied as a function of the Bjorken variable x = Q 2 /2Mu and the squared 
four- momentum transfer Q 2 = —q 2 , where M is the nucleon mass and v is the pho- 
ton energy in the laboratory frame. Accurate experimental data are now available 
for a number of nuclear targets over a wide kinematical range, 5 • 10~ 5 < x < 0.8 and 
0.03 GeV 2 < Q 2 < 200 GeV 2 . The data show nontrivial nuclear effects over the whole range 
of Bjorken x. At x < 0.1 one observes shadowing, i.e. a systematic reduction of the nuclear 
structure function F^ with respect to A times the free nucleon structure function . A 
small enhancement of the ratio R = FA j AF^ is seen at x ~ 0.2 and a pronounced dip 
occurs in that ratio at x ~ 0.5. Finally for x > 0.7 a large enhancement of R is observed. 

Numerous models have been proposed to explain these effects (for recent reviews see e.g. 
refs. 0-10]). So far, most theoretical models for nuclear DIS give separate descriptions of 



the regions of small x < 0.1 and large x > 0.2. The physical reason for such a division 
becomes apparent in the space-time analysis of the DIS process. In the laboratory frame 
the interaction of the virtual photon with the target can proceed in two possible ways (see 
Fig-©: 

(I) the photon is absorbed by a quark or antiquark in the target which picks up the large 
energy and momentum transfer; 

(II) the photon converts into a quark-antiquark pair which subsequently interacts with the 
target. 

An analysis of the contributions (I) and (II) reveals that the second mechanism dominates 



at small x <C 0.1 (see e.g. the discussion in [pl] -[T^|). At 0.1 < x < 0.3 both processes 
(I) and (II) are important, whereas the region x > 0.3 is governed by mechanism (I). In 
order to understand the implications of this for nuclear targets let us consider characteristic 
space-time scales. Mechanism (I) has a characteristic scale which is determined by the size 
of the nucleon and does not depend on x. For mechanism (II) the propagation length A of 
the qq (or hadronic) fluctuations of the photon in the laboratory frame is A ~ (Mi)" 1 . For 
x < 0.1 this propagation length becomes larger than the average nucleon-nucleon distance in 
nuclei. As a consequence deep-inelastic scattering from nuclear targets at small x involves 
the coherent interaction of the pair with several nucleons in the nucleus. This leads to 



nuclear shadowing. Models of the shadowing effect in the laboratory frame |l I - |13[|15| -|19 
usually include mechanism (II) only. 

At large x > 0.2 effects resulting from coherent multiple scattering in the nucleus are 
not important since the space-time scales of both mechanisms (I) and (II) are of the order 
of the nucleon size. In this region of x the virtual photon interacts incoherently with bound 
nucleons. Model descriptions of nuclear structure functions in the large x > 0.2 region 
p0^p7| (for a review see ||J2£J) usually start out from the impulse approximation in which 
both mechanisms (I) and (II) are taken into account by using a phenomenological nucleon 
structure function. Results of such calculations show that nuclear binding and Fermi motion 
are responsible for the observed "old" EMC-effect at large x. 
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The purpose of the present paper is to develop an approach to DIS which is based on 
a unified description of processes (I) and (II). In Section |TJ we develop a relativistically 
covariant formalism which incorporates the standard parton model but also permits us to 
include non-perturbative features which turn out to be important at small x. Our starting 
point is a general representation of structure functions in terms of dispersion integrals over 
quark spectral densities (Section |IIB| ). In Section |l| we develop a model for the quark 
spectral densities which separately reproduces the valence and sea quark parts of the free 
nucleon structure functions. In Section [TV] we discuss nuclear structure functions. We find 
that shadowing at small x has a scaling contribution from mechanism (II) (independent of 
Q 2 ) which turns out to be insufficient, however, to reproduce the empirical A-dependence. 
At this point our results differ from those of ref. ||16|| . We conclude that the propagation 
of strongly correlated qq pairs through the nucleus, partly in the form of vector mesons, is 
important to reproduce the observed shadowing effect (Section [IV A| ). For x > 0.2, where 
nuclear binding and Fermi motion are relevant, we discuss the limitations of the standard 
convolution model (Section [IV B| ). We point out that there is no reason to ignore, as is usually 
done, the dependence of the structure functions on the invariant mass of the nucleon, p 2 . 



II. FRAMEWORK 

According to the optical theorem inclusive inelastic scattering of an electron/muon on a 
nucleon or nucleus can be described in terms of the forward scattering of a virtual photon. 
The amplitude for forward Compton scattering is 

T^(P, «) = -</ *t t**(P\T &(e)j*(0)) \P), (1) 

where q and P are the photon and target momenta, respectively. The electromagnetic 
current is denoted by j^. In what follows we shall discuss the scattering from an unpolarized 
target and assume that the average is taken over target polarization in eq. ([!]). In this case 
there are only two independent terms in the Compton amplitude, 

T, u (P,q) = T 1 (x,Q 2 ) (^-9,u)+^^L,L v , (2) 

where Q 2 = —q 2 , = P^ — q^P-q/q 2 , and x = Q 2 /2P-q is the Bjorken scaling variable (we 
use the normalization (p\p') = (27r) 3 2p S ( - 3 \p — p') both for fermions and bosons, so that 
the amplitudes Ti and T 2 are dimensionless) . The structure functions F\ and F 2 are given 
by the imaginary parts of the scalar amplitudes in (fj): 

F 1>2 (x,Q 2 ) = -^-lmT lt2 (x,Q 2 ). (3) 

It is commonly assumed that in the region of high momentum transfer, Q 2 3> M 2 , the 
main contribution to the Compton amplitude comes from the diagram, Fig. |2|, in which the 
virtual photon couples to the quark current. Let us examine this contribution in detail. To 
leading order in Q 2 and in the axial gauge the Compton amplitude reads 
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ft — fi + ie 



7m 



(4) 



where the sum is taken over flavor and color degrees of freedom of the interacting quark 
which carries electric charge e a . For simplicity we have dropped the quark mass in eq.(|j). 
Here A a (k, P) is the Fourier transform of the correlator of the quark fields in the target, 



A a (k,P) = -ij d^e lk <(P\T(MOM0)) \P). 



(5) 



The structure functions can easily be found from eq.([|) by applying appropriate projection 
operators. Neglecting terms of order 1/Q 2 we find for the structure function F 2 , 



F 2 (x) 



X 



J2 e l (?a0) +q a {x)), 



where 



and 



fa(x) 



Qa{x) = f a (x), 
Qa(x) = -fa(-x), 

d 4 k Tr(^A a (fc,P)) 



(27T) 



2P-q 



5 x 



k-q^ 

~P~q, 



(6) 

(7a) 
(7b) 

(8) 



The two terms in eq.(^) correspond to the direct and crossed terms of the Compton ampli- 
tude. It follows from eq.(f|) that the structure function F\ is not independent; it is given by 
the Callan-Gross relation, Fi{x) = 2xFi(x). 



A. Parton model 

Equation (§) gives a Lorentz covariant representation of the quark distribution function 
which can be used in any reference frame. Here we first demonstrate that eqs.(||§) recover 
the familiar result of the parton model. Let us choose a reference frame in which the target 
moves with a large momentum |P| — > oo. In this frame the function q a (x) can be identified 
with the momentum distribution of quarks with flavor a in the target, and q a {x) is the 
corresponding antiquark distribution. In order to see this we introduce a coordinate system 
such that the momentum transfer is q = (0, 0±,Q) (with Q = y/Q 2 ). The hadron moves 
with momentum P 3 = —Q/2x in the direction opposite to the three-momentum transfer q. 
Then eq.(§) becomes 

f , , ■ f d 3 k ( h\ fdk Tr( l3 A(k,P)) 

where we have suppressed quark flavor and color indices for simplicity. It can be shown that 
73/-P3 = lo/Po i n eq-©) U P to terms of order M 2 x 2 /Q 2 . We perform the /^-integration by 
closing the integration contour in the upper half-plane and obtain: 



-i J ^ Tr(7 ° 2 A ^ P)) =fd^r'e-^moy)mr)). (10) 



Here we have used translational invariance and introduced an additional space integration. 
The brackets in eq.(0) denote (■■■) = (P\ ■■ ■ \P)/(P\P). The r.h.s. of eq.fllO]) is in fact the 
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momentum space density of quarks in the target. In order to clarify this connection further 
we expand the quark fields in terms of plane wave spinors, 



t cfi h 1 

r ) = E / Jo~YA~i == °) U ^ ^K kr +^, °)v(k, a)e" 4kr ) , (11) 



where u(k, a) is the Dirac spinor of a quark with momentum k and polarization a, and 
v(k, a) is the corresponding antiquark spinor. The operators a(k, a) and b(k, a) acting on 
any physical state probe the momentum distribution of quarks and antiquarks in that state: 

N(k) = ^(a)(k, a)a(k, a)), (12a) 
N(k)=J2(b\k,a)b(k,a)). (12b) 

Using the commutation relations between a, and b, W and the orthogonality properties 
of spinors u(k, a) and v(k, a), we find from eqs. (p|,|T0D , 



f(x) = q(x)-q(-x), (13) 

where q(x) and q(x) are the quark and antiquark distributions as functions of the target 
momentum fraction, 

?(*)= f^fx-^) N(M, (14a) 



(2tt) 3 V Pi 



3 , 



^=iwfA x -^) m - (i4b) 

The functions (0) have the usual simple interpretation: they represent the probability 
distributions of quarks (q(x)) or antiquarks (q(x)) which carry a fraction x of the target 
longitudinal momentum. It is known from the parton model analysis that the proba- 
bility to find a parton moving backward (x < 0) vanishes as |P| — ► oo. Also, momentum 
conservation does not permit partons with x > 1. Therefore the functions g(x) and g(a;) 
vanish outside the physical interval < x < 1. 



B. Structure functions as dispersion integrals 

Our primary task is to study deep-inelastic scattering off nuclei. In this case the simplicity 
of the parton model is lost, because no reliable approach exists to deal with nuclear systems 
in the infinite momentum frame. For this purpose the preferable frame of reference is the 
laboratory system. In the present paper we describe the distribution function f(x) using 
the analytical properties of the quark correlator (§). This method preserves relativistic 
covariance and can therefore be used in any frame. 

Following ref. Q we assume that the quark correlator A is an analytic function of the 
variables s = (P — k) 2 , u = (P + k) 2 and k 2 . For real s and u the quark correlator has 
a right-hand cut in the variable s, a left-hand cut in the variable u and singularities for 
k 2 > 0. In order to make use of these analytical properties of A in the loop integral (§), it is 
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convenient to parameterize the loop momentum k in terms of the external momenta P and 
q (introducing the Sudakov variables): 



k = aP + (3q + k ± , (15) 

where q' = x P + q, and k± is a two-dimensional vector (with k\ < 0) perpendicular to both 
P and q. 

The integration with respect to (3 can be done using the analytic properties, just men- 
tioned, of the quark correlator. As a result one finds that the distribution function f(x) 
vanishes outside the physical interval \x\ < 1 as it should. For < x < 1 the distribution 
function f(x) (or q(x)) is given by a dispersion integral in the variable s along the right- 
hand cut. For — 1 < x < the w-channel cut is relevant. In order to represent the structure 
functions in this way we introduce functions pn t L for the imaginary parts taken along the 
right-hand cut (R) and left-hand cut (L) as follows: 

2 Im^Tr^A^P)) 

p R {s, k ,a) = — — , (16a) 

Z7T P-q 

2 Im L Tr(^A(fc,P)) 

p L (u, k , a) = — — . 16b 

Z7T P-q 

In terms of these functions the quark distribution q(x) and the antiquark distribution q(x) 



arc 



g( z ) = T~— I ^ftlt PR{s,k 2 ,x), (17a) 



1-xJ 2(2?r 



Q( x ) = T— I k t Pl{u, k 2 , -x), (17b) 



1-xJ 2(2tt 

where we have suppressed the flavor indices again. The squared quark four-momentum k 2 



is 



k 2 = x f-i- + M 2 ) + - *L (18) 
Vx — 1 / 1 — x- 

in eq.(^a), and an analogous expression holds with s replaced by u for the antiquark 
distribution q(x) in eq.(|T7|b). 

The spectral densities ([TB]) can be written in terms of spectral sums over a complete set 
of intermediate states inserted between the two quark field operators in eq.(H). In order to 
write the spectral representation in a more explicit form we introduce matrix elements of 
the quark field operator taken between the nucleon and some intermediate state, ^ n (K) = 
(K.,n\ip(0)\P) and ip m (K.) = (P\ip(0)\m, K). Intermediate states are labeled by their total 
momentum K and other quantum numbers denoted by n or m. In terms of the amplitudes 
ip n and ipm the spectral densities (|1^) can then be expressed as follows: 

Pn(s,k 2 ,a) = -i-J^JP-kMVnCP-kWa " M ')' ( 19 ) 
- p L (u, k 2 , -a) = 53? m (P-k)^ m (P-k)tf(« - Ml). (20) 
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Here M 2 and are the invariant masses of the intermediate states. Note that eqs.(|19"D 
and (|20f) correspond to two different time orderings of the quark operators in (||). 

The basic assumption is now that the spectral densities vanish at large k 2 so that in- 
tegrals in eqs. (|T7|) are convergent and dominated by the region of k 2 ~ m 2 , where m is a 
characteristic hadronic mass scale. It follows from eq.(^) that the behavior of the quark 
distributions at x — > 1 is given by the asymptotics of the spectral densities at k 2 — > — oo. At 
small x the k 2 is finite even for large s ~ m 2 /x. Therefore for small x the integral in eqs.flTTD 
is sensitive to the high-energy parts of the spectral densities. In the region of intermediate 
or large x > 0.2 the region of finite s ~ rh 2 is of major importance. 

The spectral representation ( |I9| , pO"D offers a convenient way to separate contributions 
from mechanisms (I) and (II). To see this, consider the amplitude ip n in the laboratory 
frame, with target momentum P = (M, 0). The quark field operator acting on the target 
state can either annihilate a quark in the target or create an antiquark. In the former case 
the amplitude ip n oc (— k, n\a{k)\P) describes the absorption of a virtual photon by quarks 
with momentum k. This contribution corresponds to the mechanism (I). The contribution 
from the antiquark part of the ip-o^etaXoT, oc (— k, n\tf(— k)\P), corresponds to the 
mechanism (II). This part of the amplitude ip n describes the "external" antiquarks from qq 
fluctuations of the virtual photon and their interaction with the target. The amplitude ip m 
also has two parts, one associated with contributions from antiquarks bound to the target 
(V^ oc (— k, m\b(k)\P}) and the other one from "external" quarks coming from the photon 
wave function, (^ oc (— k, m\a\— k)\P)). 



III. QUARK SPECTRAL DENSITIES AND NUCLEON STRUCTURE 

FUNCTIONS 

In this section we construct nucleon structure functions using a model for the quark 
spectral densities pi and pr. We further elaborate the concept, appropriate in the laboratory 
frame, that the full quark spectral function can be divided into two parts corresponding to 
mechanisms (I) and (II), as illustrated in Fig. |l| 

The part of the spectral density that describes mechanism (I) is proportional to the prob- 
ability to find a quark with four-momentum k in the nucleon. The characteristic momenta 
of quarks bound in the nucleon are of the order of the nucleon mass, k ~ M. Therefore 
the main contributions to the spectral densities from mechanism (I) come from the region 
\k 2 \ ~ M 2 and s ~ M 2 . This kinematical region determines the behavior of structure 
functions at x > 0.2. 

On the other hand, contributions to the quark spectral densities from mechanism (II) 
rise with s (or u). In fact, the matrix element (k, n\b\k)\P) describes scattering of an 
antiquark with momentum k from the nucleon, with transition of the system to the final 
state |n, k). The possible antiquark momenta k are determined by the wave function of the 
photon and can be as large as the photon momentum q. Therefore, the invariant mass of the 
antiquark-nucleon system is large, s ^> M 2 . Due to unitarity, the sum over all states n in 
eq. (|l9D will be proportional to the antiquark-nucleon forward elastic scattering amplitude. 
It is known from Regge theory that imaginary parts of amplitudes rise with energy as s ap 
where ap is the intercept of the Pomeron. Therefore contributions to the spectral density 



S 



from mechanism (II) grow with energy and dominate at large s. This mechanism determines 
the small x part of structure functions. 



A. Model for quark spectral densities 

With this discussion in mind we now develop the following simple model for the quark 
spectral densities. We introduce a parameter so which separates the full spectrum into a 
low-energy (s < sq) part and a high-energy (s > sq) part. We assume furthermore that the 
low-energy part of the spectrum is dominated by mechanism (I), while the high-energy part 
is given by mechanism (II) (for illustration see Fig. [3]): 

p(s, k 2 ) = p®(s, k 2 )d(s -s) + pM(s, k 2 )6(s - s ). (21) 

For the low-energy part of the spectrum (s < Sq) we make the following ansatz: 

$(s,k 2 ) = <S>(k 2 )5(s-s), (22) 

together with pi = in this region, which implies that we neglect contributions to spectral 
densities coming from antiquarks in the nucleon. This choice can be motivated within a 
constituent quark picture of the nucleon. In this case ~ |M is an average mass of 
the residual two-quark intermediate state, and <£>(/c 2 ) is the squared momentum space wave 
function of the constituent quark. Here we will not confine ourselves to some particular 
model but rather choose the spectral density at small values of s in such a way that we 
reproduce the measured valence quark distribution. The contribution from eq. (|22"D to the 
deep-inelastic structure function is, 

g (I) W = T ^/_J dk 2 $(k 2 ), (23) 

where ^(s, x) is the maximum value of the squared quark four- momentum for given values 
of s and x, 

k 2 m ^s,x)=x(-^ + M 2 y (24) 

In order to describe the region of large s > Sq we introduce the quark-nucleon forward 
scattering amplitude T{k, P) as follows: 

A(k,P) = - -T(k,P)- -. (25) 

The quantity relevant for the calculation of the quark spectral densities ( |16D is Tr(7 M A). In 
terms of the amplitude T this trace can be written as follows, 

Tr( 7M A(A;,P)) = (k 2 - m\y 2 [Ak^T + [m\ - k 2 ) Tr( 7/J T)] , (26) 

where T = ^Tr[(^ + m q )T(k, P)] is the amplitude averaged over the quark spin. When 
examining the Dirac structure of T, we find that only a scalar term and a term proportional 
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to the 7-matrix contribute to (p6|). One can neglect the scalar term because its contribution 
to the amplitude T is proportional to the quark mass m q . The leading contribution to the 
amplitude T, the one which rises with s, comes from the term proportional to 7 -P. Based 
on these arguments we write the amplitude T as 

T(k,P) = C(s,u,k 2 )f, (27) 

where C is a Lorentz invariant function of s, u and k 2 which is related to the spin-averaged 
amplitude T as follows: 

T(s,u,k 2 ) = ^(u- s)C(s,u,k 2 ). (28) 

We emphasize here that the amplitude T(s, u, k 2 ) describes both the qN- and the qN- 
scattering channels. In the s-channel T(s, u, k 2 ) coincides with the antiquark-nucleon scat- 
tering amplitude TqN(s,k 2 ), while in the w-channel T(s,u,k 2 ) gives the quark-nucleon am- 
plitude T q n(u, k 2 ). 

We are now prepared to calculate the contribution from mechanism (II) to the quark spec- 
tral densities flTED and to the quark and antiquark distributions fllTf) . Using eqs. ( p6] , p7| , p8D 
we obtain :[] 

rrrw s 1 f°° , f k ^ s ' x ^ 9 lmT m (s, k 2 ) ( m 2 -k 2 \ f 

^ x)= wL ds L dk w^{- x+ 7^^)> (29a) 

-fii)/ x 1 f°°, fkl,Au, x ) lmT gN (u,k 2 ) ( m 2 -k 2 \ 

where we have introduced an integration over the squared quark four-momentum k 2 instead 
of integration over transverse momentum k±. 

In what follows we shall consider the structure function F 2 of an isoscalar nucleon, 
F2 = |(-F|2 + F-i)- Isospin symmetry implies that F^ 1 is proportional to the flavor singlet 
combination of quark and antiquark distributions, 

F 2 N (x) = ^-xj:(qa(x) + q a (x)). (30) 

(Here we have assumed that the difference between strange and charmed sea is negligible 
small). In our model the quark and antiquark distributions are given by 

E Mx) + q a (x)) = q®(x) + NfN c (q^ l \x) + q^(x)) , (31) 

a 

where q^ is given by eq.(^3|) and represents the sum of quark distributions of different 
flavors due to the mechanism (I). The quantities g^ 11 ^ and q^ are the quark and antiquark 



1 As compared to the spinless case discussed in [p|, the spin 1/2 distributions have a generic 
factor of 2 which reflects the number of spin degrees of freedom, and the term in brackets under 
the integral replaced the factor x. 
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distributions related to mechanism (II), averaged over flavor and color. They are given by 
eqs.( p9|) , where T qN and T qN are quark-nucleon amplitudes averaged over quark spin, flavor 
and color. We have Nf = 4 and N c = 3. 

The valence quark distribution is measured in neutrino scattering in terms of the struc- 
ture function F^(x). In our model it is given by 

F^(x) = q m (x) + NjN c (x) - (x)) . (32) 
This structure function is normalized to the number of valence quarks in the nucleon. 



B. Quark-nucleon amplitude 

In order to specify and we note that the amplitudes T q ^ and T q ^ can be 
connected to observable proton-proton and antiproton-proton forward scattering amplitudes. 
We recall the well known phenomenological fact that total hadronic cross-sections at high 
energies are proportional to the number of constituent quarks in hadrons |[29|| . Hence the 



forward proton-proton amplitude can be written in terms of the quark-proton amplitude in 
the laboratory frame as: 

T„(S) _ 2 ( T qp (yS,k*) \ (33) 

Here S is the squared proton-proton center of mass energy, k is the four-momentum of a 
constituent quark in the target and y is the fraction of the target light-cone momentum 
carried by the constituent quark. A similar equation relates the antiquark-proton amplitude 
Tq P to the antiproton-proton forward scattering amplitude Tp p . In the laboratory frame 



y = (ko + ks)/M, k$ = M — \jm 2 R + k 2 , where mp is the mass of the spectator system. The 
averaging in eq. (|3"3]) is performed over the spectral function of the constituent quarks in the 



proton target. In order to estimate typical values of y and k 2 we assume mp « |M. This 
gives average values y ~ 0.3 and k 2 = (k^ — k 2 ) « —0.1 GeV 2 . 



The constituent quark-nucleon scattering amplitudes in (|33| ) might in principle be dif- 
ferent from the quark-nucleon amplitudes entering in eqs. fl2"Dp . However, at small val- 
ues of x, the quark-nucleon center of mass energy s in eq.(|29|) will be of the order of 
s ~ lGeV 2 /x. As a consequence the typical formation time of a constituent quark will 
be Tp ~ \k \/m 2 ~ s/2m 2 M ~ 1/Mx. This is comparable to the propagation length 
A ~ 1/Mx of the qq pair in the photon wave function. We can therefore assume that at 
small x the quark-nucleon amplitudes which enter in the quark distribution functions can be 
approximated by the constituent quark-nucleon amplitudes determined by hadron-nucleon 
scattering. 

Through eq.(^) the s-dependence of the quark-nucleon and antiquark-nucleon ampli- 
tudes T qN and T qN is fixed. Above the resonance region the pp- and pp-amplitudes are well 
reproduced by Pomeron exchange with an intercept ap = 1+e, and by exchange of two Regge 
trajectories corresponding to the p- and the a 2 -mesons with intercepts a p = a a2 = ap ~ 1/2. 



The forward scattering amplitudes can then be written as [ pT|j : 
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T PP (S) = R P S^(i + tan ™) + S°« (iR A - ^ + ^ cosna R \ 

2 V sni7ra# / 

T PP (S) = R P S^(i + tan + ^ - ^ + ^ 008 , ( 3 4b) 

were Rp is the residue of the Pomeron while R^ and R A stand for the sum and the difference 
of residues of the p- and ct2-trajectories. In the following we shall use the best fit parameters 
of ref. H: 



e = 0.0808, R P = 21.70 mb/GeV 2e , 
a R = 0.5475, R A = 56.08 mb-GeV (1 ~ afl) , 
Rx = 98.39 mb-GeV (1 - afl) . 

The amplitudes T q ^ and TgN should depend not only on s, but also on the squared quark 
four- momentum k 2 . We assume that the fc 2 -behavior of the Pomeron and Reggeon parts of 
the amplitudes is given as: 

T qN (s } k 2 ) = g P (k 2 )T q p N (s } 0) + g R (k 2 )T q R N (s, 0), (35) 

with functions gp,n{k 2 ) for which we chose the following ansatz: 

g PtR {k 2 ) = (1 - k 2 /A 2 PtR r n ^, (36) 

with momentum space cutoffs A P , A R and exponents np and n R . 

C. Nucleon structure functions at large Q 2 

In this section we determine gp{k 2 ) and g R (k 2 ) and other remaining parameters such 
that we reproduce the measured structure functions in the scaling region. This involves the 
antiquark distributions q(x) (eq.(^) in our model) together with the structure functions 
F 2 N (x) andF 3 N (x) (eqs.®©). 

The valence quark distribution is mainly given in terms of the momentum distribution 
$(A; 2 ) in eq.(|22|) which still needs to be specified. We use the ansatz 

§(k 2 ) = (l-k 2 /A 2 v )- nv J (37) 

with a suitable cutoff parameter Ay and a characteristic exponent ny. 

From eqs.(p"7|,18) it follows that the asymptotic behavior of the quark spectral functions 
at k 2 — > — oo determines the structure functions at x — > 1. Indeed, our model gives 

q®(x -> 1) oc (1 -x) nv_1 , (38a) 
qW{x -»• 1) oc (1 - x) np ' R+1 . (38b) 



The quark counting rule [ 29 , 33 , 34 1 for the valence distribution requires q^(x — > 1) oc (1 — x) 3 
which fixes ny = 4. Furthermore, the sea quark distribution should approach zero at x — * 1 
with a high power in (1—x). We find that the Pomeron and Reggeon exponents np = n R = 4, 
which correspond to a (1 — x) 5 behavior at large x, give a good fit to q(x). 
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The Regge parameters of the quark-nucleon and antiquark-nucleon amplitude in eqs.(p9[) 
are already fixed at an averaged squared momentum \k 2 \ ~ 0.1 GeV 2 by employing eq.(^). 
Note that the overall magnitudes of the type (II) quark and antiquark distributions ( p9[ ) are 
set by the qN and qN cross sections. Using eq.(^) we find ~ cr g N ~ \ a w ~ 13 mb at 
high energy. The remaining scales enter through the cutoffs in gp,B,(k 2 ) and Q(k 2 ), together 
with s and sq which separates high and low energy parts of the quark spectral density 
These parameters are determined by fits to the nucleon structure functions at large 
Q 2 . In our analysis we use recent NMC data for F 2 (x) and neutrino data [^,0 f° r 



F^(x), xF^(x) and q(x) at average momentum transfer Q 2 « 5 — 10 GeV 2 . Our result 
shown in Fig. uses 



A 2 V = 1.2 GeV 2 , 

A| = 2.5GeV 2 , (39) 
= 4.0 GeV 2 , 

together with s = 2s = 4 GeV 2 . (40) 

The global view in Fig. |] is supplemented by more detailed comparisons with data at 
large x (Fig. ||a) and at small x (Fig. ||b). In particular, Fig. |]a demonstrates that the 
exponents ny = rip = n R = 4 are properly chosen to reproduce tails of the distribution 
functions, while Fig. [|b focuses on the behavior at x < 0.1 where sea quarks begin to 
dominate F 2 . 

A remark concerning the value of s is in order. This parameter roughly corresponds to 
the average squared mass of the residual quark-gluon system, with one quark removed from 
the nucleon. The value of s depends on the scale at which we study the system. For example, 
in the bag model (or in a constituent quark model) i/Ib^ ~ §M, where M is the nucleon 
mass. This corresponds to the quark distributions at a low resolution scale Q 2 < 1 GeV 2 , 
which then has to be evolved to the momentum transfer Q 2 at which the experimental data 
are taken (e.g. j38| ). In terms of the quark spectral densities ( |19"1J2"0|) the Q 2 -evolution effect 



modifies the spectrum of intermediate states by adding radiative corrections which generate 
gluons and qq pairs. This shifts the quark spectra to larger values of s. In our approach we fit 
to the measured structure functions at large Q 2 and effectively incorporate the Q 2 -evolution 
effect in the parameter s, averaged over the Q 2 range of the experimental data. 



D. Structure function at small Q 2 and small x 

So far we have only discussed the scaling region, Q 2 > 5 GeV 2 . New phenomena enter at 
small Q 2 and small x. This is also where shadowing effects show up in the nuclear structure 
functions F^, so that special attention is assigned to this region. 

At small Q 2 < 1 GeV 2 the scaling behavior is violated. Furthermore conservation of the 
electromagnetic current requires that F 2 must vanish as Q 2 goes to zero. We therefore need 
a model which provides a smooth transition from the scaling to non-scaling regions. Here 
we discuss a model based on the generalized vector meson dominance (GVMD) ideas ||TT 



In the GVMD approach the structure function F 2 is expressed in terms of a dimensionless 
spectral function n(/i 2 ) of hadronic states which couple to the photon: 
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FfW^fVfiEEL^S). (4!) 

Here <jn{p 2 , S) is the cross section for scattering of a hadronic state with mass /J from the 
nucleon, and S = M 2 + Q 2 (l/x — 1) is the total squared center-of-mass energy of the virtual 
photon-nucleon system. The distribution II represents the spectrum of correlated quark- 
antiquark pairs and multi- meson states with spin and parity J n = l~. It separates into a 
low-mass part dominated by the vector mesons and a high-mass continuum part Il c : 

i% 2 )= E V-^ + ne^W-^), (42) 

v=p,u,<t> \9v J 

where the continuum starts at fio ^ 1 GeV, just above the meson mass. For the vector 
meson masses my, their coupling constants gy and cross sections ovn we use standard values 
summarized in Table I: 



V 


m v (MeV) 


<#/4tt 


<tvn (mb) 


p 


768.3 


2.38 


27 


UJ 


782.0 


18.4 


27 




1019.4 


13.8 
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Table I: Vector meson masses my [0, coupling constants gy |J40 
and cross sections ovn H4l|. 



The vector meson part is then 

^' VM = " E ( — ) [~r^\ -vs. (13) 




It dominates the structure function at small Q 2 < 1 GeV 2 . On the other hand, at large 
momentum transfer Q 2 3> m v the vector-meson contribution fl43| ) vanishes as ~ m v /Q 2 . In 
this region large masses fi 2 ~ Q 2 in the Il c part of the spectral function fl42|) take over. At 
large Q 2 ^> /Xq this piece leads to a structure function with proper scaling behavior. 

With this consideration in mind, we make the following ansatz in order to interpolate 
between the regions of small and large Q 2 : f\ 

F?(x, Q 2 ) = F^ix, Q 2 ) + F ^ {x) . (44) 

For the asymptotic part F^' as (x) of the structure function we use our model as described 
in the previous section. The parameter Ql is expected to be of order /j 2 , ~ 1 GeV 2 . Scaling 
occurs for Q 2 > Q^. 



Eq.(44) is similar to a structure function model discussed in |42] 
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In Fig. |6| we show the characteristic behavior of F 2 N (x, Q 2 ) at a small value of the Bjorken 
variable (x = 0.01) as a function of Q 2 . We find that Qq = 2 GeV 2 gives a good description of 
the data. One observes that vector mesons dominate at Q 2 <C 1 GeV 2 , whereas the scattering 
of uncorrelated qq pairs governs F 2 for Q 2 > 2 GeV 2 . At Q 2 = 1 GeV 2 , roughly one third 
of F 2 at x — 0.01 is due to vector mesons; at Q 2 = 5 GeV 2 they still contribute about 10%. 
This is in agreement with results obtained within the framework of the generalized vector 
meson dominance model UTTf . 



IV. NUCLEAR STRUCTURE FUNCTIONS 

We now turn to our main theme, namely DIS on nuclear targets and the structure 
function F 2 (x,Q 2 ). 

Let us first discuss the region of small x. We have pointed out that deep inelastic 
scattering at small x (as seen from the target rest frame) proceeds dominantly via the 
mechanism (II), and we now investigate nuclear effects based on this observation. As we 
already mentioned, for x < 0.1 the propagation length A ~ (Mi) -1 of qq- fluctuations in the 
photon wave function exceeds the average distance d ~ 1.8 fm between bound nucleons in 
the nucleus. The propagating qq pair can then interact coherently with several nucleons. 
This multiple scattering effect becomes significant when A is larger than the quark mean 
free path I = (pcjqAr) -1 , where p is the average nuclear density. In terms of the x variable 
the last condition reads x < pcr q N /M. The nuclear effects due to multiple scattering of the 
gg-pair will saturate for x < (MRa)' 1 , i.e. when the propagation length A exceeds the 
nuclear radius Ra- In this case the virtual photon converts into a qq pair already outside 
the nucleus. This pair interacts with nucleons at the nuclear surface which absorbs part of 
the incoming flux and thereby screens the inner nucleons. This is the shadowing effect which 
we study in detail in the next section. 

On the other hand coherent multiple scattering effects are not important at large x where 
A < d and the qq fluctuation has no time to scatter more then once. In this region the DIS 
process takes place mainly on a single nucleon in the nucleus. Effects due to nuclear binding 
and Fermi- motion are now important. We discuss these in detail in Section [IV B. 



A. The small x region: Shadowing 

As in our previous discussion of the free nucleon structure function F^ , we use the 
following ansatz for the nuclear structure function: 

F 2 A (x, Q 2 ) = F 2 ^ u (x, Q 2 ) + 7 ^- ?V2 F 2 A - S (x), (45) 

V t VO 

and discuss the asymptotic (scaling) and vector meson contributions to F 2 separately. 

Consider first the scaling part F 2 ,as {x). It is formally obtained from F 2 [x] in our 
model through the replacement of the quark-nucleon and antiquark-nucleon amplitudes T q N 
and TqN in eqs.(p9|) by corresponding nuclear amplitudes, T q A and T q A- For the type (I) 
distributions in eq.([H|) we simply use q^(x) = Aq$(x); the validity of this approximation 
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will be discussed in the next section. We use Glauber multiple scattering theory [[43| to 
express the nuclear amplitudes Ta in terms of the T N : 



T A (s,t) = -2is J d 2 be^'- h (A\exp N^X^-b^j - 1\A). 



(46) 



Here the averaging is performed over the nuclear wave function and the sum runs over 
all bound nucleons located at positions r.,- = (hj,Zj) in the nuclear cm. frame; q' is the 
momentum transfer, t = — q' 2 . The eikonal phase x(b) is related to the nucleon amplitude 
T N as follows: 

T N (s, t) = -2is J d 2 b e iq ' b (exp(ix(b)) - 1) . (47) 
The Glauber analysis of high-energy hadron scattering from nuclei indicates that the result 



for Ta is not sensitive to iVW-correlations and other details of the nuclear wave function ]44 
We can therefore evaluate the nuclear matrix element in eq. ([4"6"D in a simple approximation 
assuming that the squared nuclear wave function is given by the product of Gaussian one- 
particle densities: 

|vJ/ A ( ri ,...,r A )| 2 = nM 2 )- 3/2 exp(-r|/it: 2 ). (48) 

The parameter R 2 = \R\ is fixed by the nuclear root mean square radius R A = 1.12 A 1 ' 3 fin 
as determined by electron scattering data. Using this wave function one can easily calculate 
Ta in terms of T/v, with the result: 



T A (s,k 2 ) = T N (s,k 2 )Y,^r RSr^ 2 • (49) 




The effective number n e s of terms which contribute to the sum (|49|) can be estimated as 
the average number of rescatterings of a classical particle moving along the diameter of the 
nucleus: n e g = 2R/1 where I is the quark mean free path, and I ~ 3 fm for a quark-nucleon 
cross section a q N ~ 13 mb. Therefore the triple scattering term j = 3 practically saturates 
the multiple scattering series for nuclei up to A ~ 100. In the actual calculations we keep 
four terms in eq. (|4"9"D . 

The calculated ratios F^(x)/AF^ (x) are shown as dashed curves in Fig. |^ for several 
nuclei. Comparing these results with the recent NMC data 0, we find that roughly only half 
of the measured shadowing effect can be explained in this way. In other words, mechanism 
(II) alone, with quark- and antiquark-nucleon interactions constrained by high energy pp 
and pp data, cannot account for all of the observed nuclear shadowing. 

We find this not surprising, for the following reason. The experimental data at small x 
are taken at relatively small values of Q 2 . For example at x ~ 0.01 the average momentum 
transfer is Q 2 ~ 1.6 GeV 2 . But at these low Q 2 it is not justified to consider only leading 
twist contributions to the structure function. The multiple scattering of strongly correlated 
qq pairs on the nuclear target now becomes important. This brings in the vector meson 
contribution p A > VM to the nuclear structure function (45|). Its form is analogous to that of 
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F 2 ,VM for the free nucleon (see eq.(|43D), where now the vector meson- nucleoli cross sections 
o"vn are replaced by corresponding nuclear cross sections ova- The latter is related to 
the former via the Glauber multiple scattering series (f49|). Together with the scaling part 
F 2 k ' as (x), which includes the nuclear interaction of uncorrelated gg-pairs, we then calculate 
the full nuclear structure function F^(x,Q 2 ) according to eq. (|45|) . 

In Fig. |7| we compare our results for R(x, Q 2 ) = F 2 A (x, Q 2 )/ AF 2 N (x, Q 2 ), including vector 
mesons, with the data of the NMC collaboration [fj]. For every x-bin, R(x, Q 2 ) is calculated 
using the corresponding average (Q 2 ) given in [fJJ. We see that the vector mesons can account 
quite well for the missing part of the shadowing. 

In summary we find that the scaling contribution to the nucleon structure function alone 
can account for only about half of the measured shadowing effect. The other half results 
from the interactions of strongly correlated qq pairs, i.e. vector mesons, with the target. 
The fact that the observed shadowing is only weakly Q 2 dependent || has now a plausible 
explanation: although the vector meson contributions vanish at large Q 2 there is still a 
sizable shadowing effect due to the interaction of uncorrelated quarks or antiquarks with the 
nuclear target. 



B. The large x region: Convolution model 

In the region x > 0.2 nuclear structure functions are commonly described within the 
impulse approximation (see Fig. |8|) ignoring final state interaction of the nucleon debris 
with remaining nuclear system.^ The Feynman diagram in Fig. || is usually written in the 
form of a convolution (see e.g. 



F 2 A (x) = fdy D N/A (y) F 2 N (x/y). (50) 

Here the structure function F 2 (x) is folded with the (light-cone) momentum distribution 
F>^/\(y) of nucleons in the nucleus. Eq.(|(]) has frequently been applied in calculations of 
Fermi-motion and binding corrections fl20|-p6|j. The convolution formalism is also used to 
evaluate meson cloud effects [47|, exchange currents corrections |23|j48|| etc. (for a review see 
ref. @). 

However, as pointed out in a recent analysis [49], a derivation of the convolution model 



(fjOp even within the impulse approximation implies several assumptions which are generally 
not justified. In this section we re-examine the derivation of eq.(50) on the basis of our co- 



variant approach developed in Section 0. It will turn out that off-shell effects are important, 
so that the simple convolution model is generally not a good approximation for F 2 . 

We start from eq.(B) which gives the nuclear light-cone distribution function Ja^a) in 
terms of the quark correlator A^{k, Pa) in the nucleus. Consider now the diagram Fig. |8|. 
It can be written as a convolution of the quark correlator Ajv(fc,p) in the bound nucleon 
and the nucleon propagator G(p, Pa) in the nucleus: 



5 An attempt to estimate the influence of final state interaction on leading twist nuclear structure 



functions was made in ref. [45,4 
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A A (k,P A ) = -i J-^^N (A N (k,p)G{p,P A )) . (51) 

Here and in the following the "hat" on Ajy and other quantities indicates their matrix 
structure in nucleon Dirac space, and the trace Tin is taken with respect to nucleon variables. 
The four-momenta P A , p and k refer to the nucleus, the bound nucleon and the quark in 
the nucleon, respectively. The quark propagator in the on-mass-shell nucleon ([5]) averaged 
over nucleon polarizations is related to An as follows: 

A N (k,p) = ~Tr N (& N (k,p)W + M)) . (52) 

Substituting eq.(|5T|) into eq.(||) we obtain for the nuclear light-cone distribution, 

d 4 p 



/aM = ~i j j^i^N (f N (x',p 2 )G(p,P A )) , (53) 
with x A = Q 2 /2P A -q and 



Here the trace Ttq is taken with respect to quark variables and x' = Q 2 /2p-q is the Bjorken 
variable of an off-shell nucleon with four-momentum p. 

Let us now examine the Lorentz structure of /jv in the nucleon Dirac space. In general 
it can be expanded in terms of a complete set of 16 Dirac matrices. However there are only 
four independent termsf] which can be constructed from the momenta p, q and the Dirac 
matrices: 

?n = fsl + fvl» + fr 13 Vap, (55) 

with 

U = 55, (56a) 
^ = 25^+2^ (56b) 

where a a p = |[7a,7/3] and fi(x,p 2 ), i = 0,1,2,3, are dimensionless Lorentz invariant func- 
tions. The coefficients in eq.(^) are chosen in such a way that the on-mass-shell distribution 
function averaged over the nucleon spin is given by: 



4 We ignore possible terms proportional to 75 and 7^75 which do not contribute to the unpolarized 
structure functions. 
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f N (x)= lim J-Tr[(j/> + M)f N (x,p 2 )] = f + f 1 + f 2 . (57) 

We see that the tensor term in fl5"5] ) does not contribute to the unpolarized structure func- 
tions. 



Substituting (|55|) into (|53| ) one finds an equation which connects the nuclear light-cone 



distribution with the functions ff. 

r d 4 p 3 

/a(»a) = -* y ^E^feg)/,^',]? 2 ). (58) 

The functions Ci(p,q) are given by traces of the nucleon propagator G with the different 
Dirac matrices in eq. (|55|) : 

C = ^TrG(p,P A ), (59a) 

& = ^Tt(G(p,P A )t), (59b) 

C 2 = -?-Tr (G(p,P A )4), (59c) 

C 3 = £^Tr(G(p,P A )a a p). (59d) 



We emphasize here that in general eq.(|58|) does not reduce to the simple convolution 



formula ( |50l) with respect to the light-cone momentum. There are two reasons for this. 
First, terms with different Lorentz structures in eq.flSop are convoluted with correspondingly 
different nuclear distribution functions. Secondly, the structure functions in the off-shell 
region depend not only on the scaling variable x' but also on the off-shell mass p 2 of the 
bound nucleon. This will now be examined in more detail. 

In the following we shall consider the nuclear structure functions as functions of the 
"nucleon" Bjorken variable x = Q 2 /2Mq instead of the "nuclear" one, x A = Q 2 /2M A q . 
The structure function F 9 function x reads: 

Ff(x) = x (F A {x) - F A (-x)) , (60) 

it f \ M t ( M \ (au 
Fa{x) = Wa U W /)- (61) 

One can easily see that the transformation (|61| ) preserves the normalization of the nuclear 
distribution function function of x: 

M A /M i 

J dxJ 7 A (x) = Jdx A f A (x A ) = 3A. (62) 

-M A /M -1 

Eq. ([58]) can be simplified if one assumes that the nucleus is a nonrelativistic system. In this 
case, as shown in the Appendix, As vanishes up to terms of order |p| 3 /M 3 if one uses the 
nonrelativistic form for the nucleon propagator. Moreover A% = for spinless nuclei. In 
the same approximation the functions Aq, Ai and A 2 are proportional to each other (see 
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eqs.(|A7]) in the Appendix). This allows us to introduce one unique nucleon distribution 
function Z?n/a which however depends also on p 2 . Using eqs. (|A7|) in the Appendix we have, 

Fa{x) = j — J dp 2 D N/A (y,p 2 )f N (x/y,p 2 ), (63) 

y>\x\ 

f N (x,p 2 ) = ^p 2 /M 2 f (x,p 2 ) + (p 2 /M 2 )h{xy) + f 2 (x,p 2 ), (64) 

^N/A(y,P 2 ) = /^W)(l + |)%-§)^ 2 -^ 2 ). (65) 

Eq . (|64"D can be identified with the light-cone distribution function of the bound nucleon. In 
eq.(|B"oD we have introduced the nuclear spectral function 

S(p) = 2tt ]T S(po ~ M - e n ) |v]/ n (p) | 2 , (66) 

n 

where the sum includes all residual nuclear states with A — 1 nucleons which carry together 
the momentum —p. All other quantum numbers are denoted as n. The nucleon separation 
energy is defined as e n = E (A) - E n (A - 1). Furthermore ^„(p) = ((A - l) n , -p|#(0)|A) 
where ^(O) is the nonrelativistic nucleon field operator at r = 0. The spectral function is 
normalized to the number of nucleons in the nucleus, 

' /V S(p) = A, (67) 



(27T) 



which guarantees the correct normalization of the nucleon distribution function in eq.(|65|). 

It is convenient to introduce also the distributions of quarks, q(x,p 2 ), and antiquarks, 
q(x,p 2 ), in the bound nucleon. These are expressed through /iv(x,p 2 ) by eqs.(^): 

Q(x,p 2 ) = f N (x,p 2 ), (68a) 
q(x,p 2 ) = -f N (-x,p 2 ) (68b) 

In terms of these distributions the structure functions of the bound nucleon are given by the 
usual parton model formula (see eq.(^])). The relation between the nuclear and the bound 
nucleon structure functions then reads: 

M A /M 

F 2 A (x)= J dyJdp 2 D N/A (y,p 2 )F 2 N (x/y;p 2 ). (69) 

X 

Let us examine this equation in more detail. The nucleon distribution (^) is strongly 
peaked around p 2 = M 2 and y = 1, with a characteristic width Ay ~ p F /M, where pp is 
the nucleon Fermi-momentum. Expanding the bound nucleon structure function in eq.(|63"D 
in a Taylor series around these points and integrating term by term, one then obtains the 
following expression for the nuclear structure function per nucleon: 

F 2 \x)/A ~ F 2 N (x) - M xF f(x) + M f F »\ x) 

+ 2 - TT P L* ■ ( 7 °) 



M V dp 2 



p2=M 2 
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Here F^'(x) and F^"(x) are derivatives of the structure function with respect to x, and (e) 
and (T) are the mean separation and kinetic energies of the bound nucleon, 

Corrections to eq.(|70|) are of higher order in (e)/M and (T)/M. One should also note that 
eq. (|70|) can safely be used for 1 — x > Pf/M ~ 0.3. In this region the condition x/y < 1 gives 
practically no restrictions on the integration over the nucleon momentum p in eqs.(|7l|,^). 
The first three terms in eq.fl70D are identical to the result of [20,23| in their discussion of 



the EMC-effect, while the last one reflects the leading contribution from the p 2 -dependence 
of the bound nucleon structure function. Let us first neglect the latter and discuss effects due 
to separation and kinetic energies. These terms are of opposite signs and the competition 
between them results in a behavior of the ratio R(x) = F^ jAF^ at x > 0.3 similar to 
that seen in the experiment. As it was first pointed out in ref. (gOj and discussed by many 
authors (see e.g. ref. || for a review), this may account for the EMC-effect at intermediate 
and large x. An important (and still open) problem in this respect is a reliable calculation 
of (e) and (T). In a simple nuclear shell model the removal energy is averaged over all 
occupied levels. One finds typical values (e) « -(20 - 25) MeV and (T) ^18-20 MeV. 
Correlations between nucleons change the simple mean field picture significantly and lead 
to high momentum (p > pp) components in the nuclear spectral function (0). This in turn 
causes an increase of the average removal energy (e). In order to demonstrate this let us 



consider the Koltun sum rule [p0| , which is exact if only two-body forces are present in the 
nuclear Hamiltonian: 

(e) + (T) = 2/i B . (73) 

Here \lb ~ —8 MeV is the nuclear binding energy per nucleon. In particular, eq.([F3|) tells 
that an increase of (T) due to high momentum components implies also an increase of 



| (e) | . We refer in this respect to a recent calculation |5T] of the spectral function of nuclear 
matter based on a variational method. This calculation shows that there is a significant 
probability to find nucleons with high momentum and large separation energies. Integration 
of the spectral function of ref. f5lj gives (T) « 38 MeV and (e) « — 70 MeV. In order to 



estimate these quantities for finite nuclei one usually assumes [25] that the high momentum 
component of the nucleon momentum distribution is about the same as in nuclear matter, 
which gives (T) w 35 MeV for a wide range of nuclei. The latter quantity together with 
eq.(|73|) leads to (e) ~ —50 MeV. It should be noted however that, even though a qualitative 
understanding of the EMC effect can be obtained using such values for (e) and (T), a 
quantitative description is still lacking. 

Let us finally discuss nuclear effects due to the off- mass shell properties of nucleons 
bound in nuclei. We note in this respect that the analysis of section |J can be applied 



also to eq.fl54) which describes the quark and the antiquark distributions of the off-shell 
nucleon. We parameterize the loop momentum in eq.fl54]) in terms of the Sudakov variables, 
k = ap + f3q' + k±, where p is the nucleon momentum in eq.([5l|) and q' = q + x'p. Assuming 
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again that the quark correlator A at is an analytic function of s = (p — k) 2 , u = (p + k) 2 
and k 2 we end up with equations similar to eqs.( |T7| , |I%D , where however the squared nucleon 
mass M 2 is now replaced by p 2 . In particular, for the quark distribution q(x',p 2 ) we have 

q(x',p 2 ) = — ^ / ds / ma *dk 2 p R (s,k 2 ,x',p 2 ), (74) 

1071" J J-OD 

where £^ ax = x' (s/(x f — l)+p 2 ) (cf. eq. ([18])) and pr is the quark spectral density given by 
an equation similar to eq. (|16D . 

We conclude from eq.(|74"D that the p 2 -dependence of the structure functions has two 
primary sources: 

1) Explicit dependence of the quark correlator Ajy in eq.(|54|) on p 2 ("dynamical" p 2 - 
dependence). This leads to a p 2 -dependence of parameters characterizing the quark 
spectral density pr. 

2) Dependence of the invariant variables s and u on p 2 ( "kinematical" p 2 -dependence) . 
This manifests itself via the p 2 dependence of k^^ in eq.(|74j). 

In this section we investigate nuclear effects at x > 0.2. Therefore we keep contributions 
to the quark spectral density pn only from mechanism (I) which dominates in the large x 
region. Recalling eq.(^) we now have 

p R (s,k 2 ,x>,p 2 ) = <S>(k 2 ,p 2 )6(s-s), (75) 

and evaluate the derivative of the bound nucleon structure function with respect to p 2 : 

dFj*(x,p 2 ) _ dqW(x,p 2 ) 
dp 2 dp 2 

^(^„y )+ />«). , 76) 

The two terms in eq. (|76|) reflect the two sources of the ^-dependence. The first one arises 
from the p 2 -dependence of fc^ ax . Its contribution to F<f in eq. ([?(]) leads to an enhancement 
of the binding effect (see Fig. ||). However if only this "kinematical" p 2 -dependence would 
be present the number of valence quarks in the nucleon would change with p 2 . Therefore 
an explicit p 2 -dependence of the quark spectral density is necessary to fix the normalization 
of the valence quark distribution. In order to evaluate this effect we assume that the k 2 - 
dependence of the function $(/c 2 , p 2 ) is the same as for the on-mass-shell nucleon, eq . (p6|) , and 
that the p 2 -dependence comes from the corresponding dependence of the cut-off parameter 
Ay, $(/c 2 ,p 2 ) = $(A; 2 ; Ay(p 2 )). We fix the p 2 -dependence of Ay in such a way that the 
number of valence quarks in the nucleon is independent of the off-shell nucleon mass p 2 , i.e. 

JLj Q 1 dxqV\x,p 2 ) = 0. (77) 

Together with eq. (ff6|) we can now calculate dAy/dp 2 . We find dAy/dp 2 ~ 0.15 at p 2 = M 2 . 
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The p 2 dependence of Ay can be viewed as an effect of changing the bound nucleon size 
Rn ~ 1/Ay in nuclei. In fact dAy/dp 2 > implies a "swelling" of the nucleon in the nuclear 
environment. For the relative change of the radius of the bound nucleon we obtain, 

SR N ISA 2 V dA 2 v M(e - T) 

R N ~ 2 A 2 V dp 2 A 2 V " 1 j 

In this way we find a very small (0.8%) increase of the size of a nucleon bound in the 
nucleus, indicating a remarkable stability of the nucleon.0 This result is compatible with 
recent findings using quite a different approach f46|| . 

In Fig. H we show a typical result for R(x) = F^(x)/(AF^(x)) calculated using eq. (ffOj) . 
The solid curve corresponds to the impulse approximation neglecting any p 2 -dependence. 
The dashed curve is the result of the full calculation including eq. fl76|) . The "kinematical" 
and the "dynamical" p 2 -dependence tend to cancel each other partly at small x < 0.3. At 
x ~ 0.5 — 0.6 the effect of the p 2 -dependence of the bound nucleon structure function is 
clearly visible and leads to an enhancement of the nuclear binding effect. 

A result (for 40 Ca) of the unified description which incorporates both the shadowing 
effect at small x and the binding, Fermi motion and off-shell corrections at large x is shown 
in Fig. [II]. While the overall pattern of the data is quite well reproduced, we see that once 
the off-shell p 2 dependence of the bound nucleon structure function is included, there is some 
room for a possible small enhancement of F^ due to nuclear pion cloud effects [ p3| , f48"| for 
0.2 < x < 0.4, which are omitted in the present calculations. 



V. SUMMARY AND CONCLUSIONS 

We have presented a unified description of deep inelastic scattering on nuclear targets 
which covers the whole region 10~ 3 < x < 1 of the Bjorken variable. Our starting point 
is a relativistic, covariant formalism which makes use of the analytic properties of quark 
correlators. In the infinite momentum frame we recover the usual parton model. In the 
laboratory system, which is the appropriate frame in which to investigate nuclear structure 
functions, this approach naturally incorporates two basic mechanisms, namely (I) scattering 
from quarks and antiquarks in the target and (II) photon conversion into a quark- ant iquark 
pair and subsequent interactions of this pair with target constituents. 

For small x, say below x < 0.1, the second one of these processes dominates and produces 
shadowing. At the momentum transfers, Q 2 , typical of the current experiments we find 
that only about half of the observed shadowing comes from the coherent interaction of 
uncorrelated qq pairs with target nucleons, the mechanism discussed in ref. \TE\. The other 



half comes from the coherent scattering of strongly correlated pairs, i.e. vector mesons. 

At larger values of the Bjorken variable (x > 0.2), scattering from the quarks of a single 
bound nucleon dominates. The leading nuclear effects are now binding and Fermi motion, 



5 In ref. |2^] the EMC effect is attributed to a 10% increase of the nucleon radius in the nucleus. 
We emphasize here that in our study of nuclear binding and off-shell effects this increase of size is 
reduced by an order of magnitude. 
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as pointed out already in previous studies. However, we find that the naive convolution 
formalism needs to be generalized to incorporate off-shell effects characteristic of bound 
nucleons. These off-shell effects are by no means small. They reduce FA(x)/A with respect 
to the free nucleon structure function F%{x) and effectively enhance the binding correction. 

We thank W. Melnitchouk and A. W. Thomas for discussions and comments. 



APPENDIX: NONRELATIVISTIC REDUCTION OF MATRIX ELEMENTS 

Consider eq . fl58|) written for the distribution function !Fa(x). For the following it is 
convenient to multiply it by x. We have: 

xF A {x) = -i /t^EAx'/,^/). (Al) 
J ( 27r ) i=0 

where A = ^ Tr G(p) (A2a) 

Ai = ^Tt{G(p)t) (A2b) 
A 2 = Tr(G(p) 7+ ) (A2c) 

As = \% Tr (p(p)°°p) ■ ( A2d ) 

Eqs. (|A2|) are written in the laboratory frame where Pa = (Ma,0), q = (qo, Ot, — |q|), 
P+ — Po + P3 an d 7+ = 7o + 73- We have used the fact that x/x' = p + /M. Also we 
have introduced G(p) = G(p, Pa)/2Ma_ where the factor IMa ensures that G is independent 
of the normalization of the nuclear state. 

Consider a non-relativistic reduction of matrix elements (A"2). We start from the nucleon 
propagator which can be written as follows: 

G(p) =-ij dte*°4(T(N(p,t)N(p,0j)), (A3) 

where N(p,t) = J d 3 r exp(— ip ■ r) N(r, t) is the nucleon field operator in a mixed (p, i)- 
representation and the brackets denote the averaging over the nuclear ground state, (••■) = 
(A\ ■ ■ ■ \A)/(A\A) (see also Section 0). 

We apply a 1/c-expansion technique to obtain an approximate solution of the Dirac 
equation for the nucleon field N in the non-relativistic limit |p|/M — > 0. Up to terms of 
order 1/c 3 the nucleon field N can be written as follows: 



*(P,*) 

TXT 



WCP,t)=C( J?S ,^ ( 1. (A4) 



We have introduced a two-component non-relativistic nucleon field The normalization 
constant C is fixed by the charge (particle number) conservation condition, 

J d 3 r A f (r)A(r) = J d 3 r ^ f (r)^(r), (A5) 



24 



which gives C — 1 — p 2 /8M 2 



Now we are prepared to calculate traces (S2). We write the nucleon four- momentum as 



p = (M + £, p). The squared four- momentum is p 2 ~ M 2 + 2M(e — T) where T = p 2 /2M is 
the nonrelativistic kinetic energy. We also introduce the nonrelativistic nucleon propagator 

G(p) = -ij dte^-*(r(*(p,t)*t(p )0 ))) ) (A6) 
Using eqs.( |TOyX6| ) we have: 



MP) = V^ (1 + M )tr ^ (p) ' (A?a) 
Mp) = j^(l + %)KG(P), (A7b) 
A(p) = (l + g)tr^(p) (A7c) 



As( P ) = i^ ] ^tr[g(p)(px<T) 3 }, (A7d) 



where the trace is taken with respect to spin. Corrections to eqs. ( |A7| ) are of order 1/c 3 . From 
the last equation we see that ^,3 appears only in the order |p| 3 /M 3 and can therefore be 
neglected at the present level of accuracy. Moreover this term vanishes identically for spinless 
nuclei, where tr^cr = 0. Substituting eqs.( [A7D into eq. ([All) an d closing the integration 
contour in the upper half of the complex p plane we arrive at eq.(p3f). 
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FIGURES 



FIG. 1. Two basic mechanisms of deep inelastic scattering. 

FIG. 2. Compton amplitude to the leading order in Q 2 . 
FIG. 3. A schematic picture of the quark spectral function. 

FIG. 4. The nucleon structure functions F^ix), F^{x) and the antiquark distribution q(x). 
Experimental data for the F^(x) (circles) are from ref. (j35|, for F^(x) (squares) and q(x) (triangles) 



from ref. [36 



FIG. 5. The nucleon structure functions F^(x), F^(x) and q(x) at large (a) and small (b) 
values of the Bjorken variable x. 

FIG. 6. The Q 2 -behavior of F^ \x,Q 2 ) at small x. The data points are from ref. [35] at 
x = 0.008 (circles) and x = 0.0175 (squares). The contribution from vector mesons (S3) is shown 
by the dashed line. The solid curve is the full result using eq.(H3). The dotted line is the result for 



F% (x = 0.01, Q ) using an empirical parameterization of the NMC data [35]. 



FIG. 7. The shadowing effect calculated for 12 C (a), 40 Ca (b) and 132 Xe (c) in comparison with 
the experimental data pL§[. The dashed curve represents a scaling part of the shadowing. The 
solid curve is the result of the full calculation including the vector mesons. 



FIG. 8. Impulse approximation for the nuclear Compton amplitude. 



FIG. 9. The ratio of the nuclear and nucleon structure functions calculated using eq.(70) with 
(e) = — 50MeV and (T) = 20MeV. The dashed curve is the result without corrections due to 
the p 2 -dependence of the bound nucleon structure function. The solid line is the result of the full 
calculation. The experimental data for 197 Au are taken from |7J. 

FIG. 10. Result of the full calculation for 40 Ca, including shadowing, binding and off-shell 
effects. 
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